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The primary goals of this thesis included developing a computer-integrated LDV 
system, utilizing the LDV system to measure the mean and fluctuating quantities of the 
flow field of an airfoil-spoiler system, and comparing the results of the experiment with 
previous hot-wire anemometry results. 

The experiment was conducted in the NASA Ames 11cm X 25 cm indraft 
tunnel. A two-component LDV system was developed and set up, and measurements 
were obtained along the upper surface of the airfoil (ahead of and behind the spoiler) 
as well as in its wake, close to the trailing edge. Mean velocity profiles and turbulent 
stresses were calculated and comparisons were made with previous hot-wire 
experiments. 

It was shown that the hot-wire technique is not as accurate as the LDV in highly 
turbulent regions of the flow field. The hot-wire results are also not reliable in the near 
wake region due to the inherent reverse flow. In other, less turbulent, regions of the 
flow field where there is no reverse flow, good agreement between previous hot-wire 


data and present experimental results was found. 
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I. INTRODUCTION 


A. SPOILERS 

A spoiler is a plate on the upper surface of a wing which can be deflected to spoil 
the flow and to cause early separation. Consequently, lift decreases and drag increases. 
Spoilers are used as effective aerodynamic control surfaces in transport aircrafts. They 
can be symmetrically deflected as speed brakes or as lift dumpers at touch down. 
When spoilers are deflected asymmetrically, they can provide substantial amounts of 
roll control in certain conditions of flight. 

The flow field of an airfoil with a deflected spoiler is complex: including flow 
separation, reattachment and vortex shedding [Ref l:p. 1]. 

The basic steady state structure of the flow field of an airfoil with a deflected 
spoiler has been measured and hypothesized by the Boeing Research Group [Ref. 2]. 
Further work into how the unsteady wake relates to the airfoil geometry and to other 
DPararmeters was undertaken by Lee, et al. In their experiments, Lee, et al., used 
hot-wire anemometers and pitot tubes to systematically measure the mean and 
fluctuating quantities on the airfoil surface, as well as in the wake, in order to better 


understand the structure of this flow field. [Ref 1] 


B. LASER DOPPLER VELOCINETER SYSTEMS 

Advanced research in experimental fluid dynamics requires a familiarity with 
sophisticated measurement techniques. Optical methods and in particular, the Laser 
Doppler Velocimeter (LDV) are now recognized as the most reliable means for 
performing measurements in complex turbulent flows [Ref. 3]. 

Since its inception in 1964, the LDV technique has been revised, developed, and 
tested by many researchers and its capabilities are reasonably well-understood. The 
method has the advantage of performing instantaneous, non-intrusive velocity 
measurements of small particles suspended in the flow. These measurements are 
relatively independent of fluid properties and the method does not require frequent 
calibration. Furthermore, its ability to measure the flow without disturbing it 1s 
especially attractive. However, because of the range of disciplines involved including 
optics, electronics, light scattering, and signal processing, the implementation of this 
method can be quite onerous. Scrious measurement errors can occur if the 


measurement techniques are not applied correctly. 


Ie 


C. THESIS GOALS 
Considering this previous work, the goals of this thesis experiment were threefold. 


First, to develop and set up a two-component LDV _ system with complete 
computer-experiment integration employed. Second, to utilize this LDV system to 
measure the mean and fluctuating velocities of an airfoil with a deflected spoiler flap, 
especially in the reverse flow region and the near wake where the fluctuations are large. 
Third, to compare the LDV data obtained with the previous hot-wire anemometry 


data. 
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Md EXPERIMENTAL ASPECTS 


A. EXPERIMENTAL SETUP 

Whe experiments were conducted in the NASA Ames I! cm X 25 cm indraft 
Wtoermmnecl rt icure 2.1 Shows the test section of this tunnel. The test section is 
equipped with slotted walls and a number of plenum chambers above and below the 
centerline. These features were used in previous experiments to give the tunnel 
adaptive walls. The adaptive wall concept was not emploved in this experiment. 
however the slotted walls did provide a means of reducing wall interference. The wind 
tunnel also ts equipped with interchangeable glass and plexiglas test section windows. 
For these experiments, the glass windows were used due to their higher refractive 
qualities which enhance the LDV signals. 

The basic model used was a Boeing advanced transport research airfoil section 
with a 6-inch chord and maximum thickness ratio of 11.3% as shown in Figure 2.2. 
This airfoil has a spoiler of 15.5% chord length hinged at 73% chord. The coordinates 
of the airfoil are given in Table I. This airfoil was selected because of the availability 
of hot-wire anemometer and pitot tube experimental data for comparison with the 


two-component LDV results. 


B. INSTRUMENTATION 
]. Pressure Instrumentation 

For this experiment, pressure information was obtained by taking advantage 
of a bank of six, 24-port Scanivalves tapped into the tunnel during a previous 
experiment. The six pressure transducers were of the differential type with a range of 
=35 psig. One side of each transducer was connected with the measured pressure 
through plastic tubing while the other side of each transducer was vented to 
atmosphere. 

In order to calibrate the transducers, one port of each Scanivalve was tapped 
to a calibration pressure generated by a small electric vacuum pump. Another port of 
each Scanivalve was tapped to the atmosphere. The calibration pressure and 
atmospheric pressure were each measured by alternate methods. These two known 
pressures and the values given for them by each transducer allowed the Scanivalves to 


be calibrated. 


ye 
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igure 2.1 NASA Ames indraft wind tunnel test section. 
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Figure 2.2 Boeing advanced transport research airfoil section. 


2. Laser Doppler Velocimeter Instrumentation 

Figure 2.3 shows the two-component, two color LDV arrangement used at 
NASA Ames. The laser 1s a four Watt Argon ion unit (Spectra Physics Model 64). 
This laser is configured with an end cavity multiline reflector which allows it to operate 
in the multiline mode. Translation stages in three axes permit positioning of the probe 
volume anywhere in the test section. 

In the multiline mode, the resultant beam is comprised of all wavelengths 
inherent to the Argon ion unit. This beam exits the laser in a fairly collimated state, 
hence an external collimator ts not required. The beam first passes through a 
half-wave plate (polarization rotator) where it becomes horizontally polarized. A set of 
Brewster angle prisms is then used to separate the primary beam into its various 
wavelengths. To obtain maximuin light transmission through the prisms, the beam 
polarization must be parallel to the plane defined by the incident and refracted beams 
[Ref. 4]. This was the reason that the initially vertically-polarized laser beam was 
passed through a polarization rotator prior to entering the dispersion prism. When the 
two prism method 1s used, the prisms also turn the beams which, with a little planning, 


may be used to advantage. At some distance from the prisms, the beams have 
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sufficient separation so mirrors can be used to direct the individual beams to the next 
optical components. The two stronger lines, green (514.5 nm) and blue (488 nm), are 
guided by means of turning mirrors to Bragg cells. The weaker, unused lines are 
“trapped” (blocked by nonreflective plates) to keep them from interfering with the 
system. Prior to reaching its Bragg cell, the green beam is passed through another 
half-wave plate in order to vertically polarize it. 

A Bragg cell is an acousto-optic device. Both the blue and green beams are 
passed through Bragg cells which split each into two equal-intensity, slightly diverging 
beams while inducing a 40 MHz modulation in one of the beams in each pair. The 
difference is that the green beam is split in the vertical plane while the blue beam 1s 


split in the horizontal plane (Figure 2.3). 
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Figure 2.3 Two-component LDV arrangement used at NASA Ames. 
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In each case (blue and green), the two diverging beams are then made parallel 
by passing each pair through its own alignment prism. From there on, turning mirrors 
guide the beams through the translation stages. They are finally focused by means of a 
condensing lens at the desired point in the test section of the wind tunnel. 

The four beams must all focus to a common point in the flow and be 
approximately coaxial. Using the steering mirrors for the appropriate beam pair, the 
focused beams are adjusted to visually overlap. Following this, various techniques 
exist for making fine adjustments to insure the beams all cross at a common point. 

An off-axis, forward scatter configuration has been selected for the receiving 
optics used at NASA Ames. This gives flexibility to the system so it can be adjusted to 
different requirements; if maximum signal strength is required, angular deviation from 
the beam axis of symmetry 1s minimized. If focusing on the smallest part of the probe 
volume is of interest, the angular deviation 1s made as large as the physical constraints 
will allow. 

The first lens in the receiving optics 1s the aperture through which the 
scattered light is collected (Figure 2.4). The second lens collimates the beam. A 
collimated beam 1s required by the polarization beamsplitter for maximum efficiency. 
The polarization beamsplitter allows the vertically polarized green beam to pass 
Straight through while the horizontally polarized blue beam is reflected upward. The 
use of polarization rather than frequency to separate the two wavelengths substantially 
increases the efficiency of the system. Each beam then passes through a line filter 
centered on the respective wavelength (blue or green) in order to reduce both the 
ambient light and the cross-talk due to impure polarization. Each line filter 1s followed 
bv a lens which serves to focus the light onto the pinhole aperture of the respective 
photomultiplier. The photomultipliers convert this light energy into electronic signals 
which are then passed on to the signal processors. 

Preamplifiers are located in the photodetector housings to further increase the 
signal amplitude and drive the signal over relatively long cables. The preamplifier 1s 
kept close to the detectors in order to increase the signal level before additional noise 
enters the signal cables. It is always good practice to amplify the signal as much as 
possible with the photomultiplier tube before going to the preamplifier which inevitably 
produces electronic noise. 

The signal processors used in this experiment were of the counter type and 


produced by Macrodyne. These processors have good resolution and superior noise 
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Figure 2.4 Computer-experiment integration. 


rejection circuitry. The Macrodyne unit is a high-speed counter with a | GlIIz clock 
designed to operate in a variety of experimental situations. It has been applied under 
single particle, multi-particle and multi-velocity conditions. In addition to the basic 5/8 
(or 10/16) signal validation criteria and oversize particle rejection, a sophisticated 
multi-level sequence validation has been incorporated. This allows the experimenter to 
select the criteria which will best yield data for a given signal-to-noise ratio, signal 
level, seeding level, Doppler frequency, and fluctuation, indicative of the test situation. 

When the particle seeding rate is low, as was the case in this experiment, the 
LDV signals become intermittent bursts. Under these conditions, the counter 
processor is the best choice as a means for obtaining the Doppler difference frequency. 
If the seeding is high, the signals will be almost continuous with random amplitude and 
phase fluctuations. The counter will still function but cannot be expected to produce 
accurate results to Detter than as 2-1 ere. 

The counter processor determines the Doppler difference frequency in the time 
domain. Timing over a selected number of cycles in the burst signal is accomplished 
with a crystal controlled “clock” oscillator. The frequency of the oscillator must be 
much higher than the Doppler difference frequency to achieve good resolution. The 
clock count is initiated when the high-pass filtered signal exceeds a set threshold level. 
A preselected number of cycles (8 or 16) are processed and used to stop the counter. 
Greater accuracy may be obtained by processing more cycles. Some processors, (€.g., 
Acrometrics) count all of the cycles in the burst and average the result over the number 
of cycles counted. This capability is especially important when the seed particles are 
polydisperse and produce a range of signal amplitudes. In such cases, there is the 
possibility of completing the preset number of cycles on the early part of the signal 
which has a lower signal-to-noise ratio and slew rate. There is also the possibility of 
reading signals from larger particles more than once. 

Although the importance of attaining good quality signals cannot be 
overemphasized, in many practical cases the signals will be noisy. Therefore, each 
burst must pass certain criteria before it is accepted. The so-called periodicity check is 
used to reject: 1) signals with insufficient cycles, 2) signals produced by two or more 
particles in the measuring volume at a time, and 3) noise. Two counters are used to 
make period measurements over a different number of cycles, usually 5 and 8. The 
ratio of the two measurements is then compared to five-eighths and if they agree to a 


specified percentage, the signal is validated. The problem with this technique is that 
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the entire signal mav be rejected because of a noisy cycle early in the burst. However, 
when the method ts coupled with the three-level validation circuitry, it can be effective 
in rejecting spurious noise. 

The three-level validation requires that a signal, Figure 2.5, exceed a postiive 
threshold level, pass through zero, and exceed a negative threshold level before the next 
zero crossing is counted as a cycle. Such logic has proven to be verv effective in 
rejecting noise. If the conditions are not met on any cycle, the counter resets and 
Starts again almost immediaiely. Thus, noisy cycles on the start of a burst need not 
cause the loss of an otherwise good signal. 

These amplified electronic signals are sent to the signal processor and 
ultimately to the data processor for conversion into velocities. 

Ideally, an LDV needs no calibration since 1t provides a direct measure of the 
particle velocity. The fringe spacing. 6, can be determined using the following formula 


where A is the laser wavelength and 6 is the angle of the intersecting beams: 


6 = A (2X sin(B 2)) 


Knowing the fmnge spacing and the time required for a given number of crossings, 
particle velocity can be determined. In practice this method can be inaccurate due io 
misalignment of the optics. The beam intersection angle. 9, is also very difficult to 
compute accurately. Therefore, it is Wise to test the entire system by measuring the 
light scattered from an object moving at a velocity that is known a priori. A rotating 
disk with a mark at a known radius ts usually the most convenient means of doing this. 
However, another method (the method used in this experiment) is to measure the free 
stream velocity using the LDV and compare this value with that computed from the 
Scanivalve data. 
3. Flow Seeding Requirements 

An inconvenient characteristic of the laser Doppler velocimeter technique ts 
the need for seed particles. The LDV actually measures the velocity of small particles 
that are assumed to move with the flow. If these particles do not have the appropriate 
size and concentration distribution, the experiment results mav be adversely affected. 
Particle sizes on the order of | micron are often quoted as adequate for mosi gas flows 


and can respond to turbulence frequencies exceeding one kHz. Smaller particles may 
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Figure 2.5 Three-level detection scheme. 
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be below the detectable limit of the instrument but will contribute to the background 
noise while larger particles that are easily detectable may not adequately track the flow. 

It was determined in this experimental situation that seeding of the flow was 
not absolutely necessary. Several runs were made with no seeding supplied. The data 
obtained in these runs was accurate as dust particles acted as the flow-following agent. 
However, it was determined that by using a flow seeding system the data rate was 
increased by a factor of four to five. The system used consisted of a paint sprayer 


spraying latex particles (on the order of 3 microns) suspended in alcohol. 


C. COMPUTER-EXPERIMENT INTEGRATION 

The experiment was fully integrated with a Data General NOVA muinicomputer 
for on-line data acquisition and analysis, Figure 2.4. A main program 
(“SPOILER.FR”) was written to. perform all tasks and record all data required during 
the experiment so that individual tunnel runs proceeded virtually ‘hands off. The main 
program called a host of other programs during each run in order to accomplish 
various tasks. These other routines are listed here: 

(1) “SCANTI.ER” - used to obtain pressure data from the Scanivalves and 
compute free stream velocity and free stream Mach number. “SCANI.FR’ 
itself calls three other routines in order to determine ambient temperature, 
calibration pressure and barometric pressure: 

(a) “TEMP.FR” - which vields ambient temperature from a Hssimenreioen ae 
via an analog-to-digital converter. The thermocouple is located 
upstream of the test section. 

(b) “CALPRES.FR” - which calculates the output pressure from the 
calibration vacuum eresiet ive pressure value 18 read by ‘tis 
program using an ana ee OS SEN converter. This value is then used 
in the calibration of the Scanivalves. 


(c) “NPARO.FR” - which yields the value for current atmospheric 
pressure from a Paroscientific pressure sensor via a digital I/O. 


(2ST VeErR =Seviich gives the current time and date in order that the data 
files being created by the main program can be labeled. 


(3) “TRELAY.FR” - to command the taskmaster to move the traverse mechanism 
to the next data collection location. 


(4) “PLOTP.FR’” - to provide on-line data analysis by allowing any of five curves 
to be plotted immediately after a given run. 
D. EXPERIMENTAL CONDITIONS 
The free stream velocity was maintained at a nominal value of 70 m/sec, 
corresponding to a Reynolds number of 7.4* 10° based on the airfoil chord. 
Experiments were carried out with the spoiler deflected 30 degrees and the angle of 


attack at 0 and 5 degrees. Locations where data was obtained are shown in [igure 2.0. 
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Figure 2.6 Data collection locations. 
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In order to determine the velocity profiles, data was acquired at varying steps 
based on the velocity gradient. Fine step sizes were taken near the airfoil for boundary 
layer profiles and near free shear layers in the wake measurements. Typically, 15 data 
points were taken in the boundary layer and 20 points in the wake. The 4-beam LDV 
arrangement precluded obtaining data closer than 0.1 inch from the airfoil surface, 
however. 

The boundary layer traverse was along the vertical axis of the test section, so the 
traverse direction is not normal to the airfoil surface as are usual boundary layer 
profiles. Angle of attack and model surface curvature make the measurement axis 


incline to the model. 


1S) 


Hl. DATA ACQUISITION 


A. FREE STREAM VELOCITY 

At the beginning of every run, the main experiment-driving program 
(“SPOILER.FR”) calls for the Scanivalve pressures, ambient temperature, calibration 
pressure and atmospheric pressure to be read. Assuming isentropic flow in the test 
section upstream of the airfoil and away from any shock waves, the free stream Mach 
number, M, can be calculated from the stagnation pressure, a and the static 


pressure, P, along the wall: 


(PPT DIY — Nac-dyl”? 


<< 
I 


For y= 1.4: 


(((P/P)2/7) = 1) x 5} 1/2 


2 
l 


Then, velocity is given by: 
V= Ma 
where 
a = /yRT 
T is obtained from the total temperature, T,,, using the following formula: 


T = Ty/{l + ((y-1/2) x M7} 
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T,, 1s measured by a thermocouple at the wind tunnel inlet. 
The above method 1s used to obtain free stream Mach number and free stream 
velocity for every run. The value for free stream velocity is used in the data reduction 


process to normalize the airfoil velocities obtained from the LDV. 


B. AIRFOIL VELOCITIES 

Processing of the Laser Doppler Velocimeter data obtained during this 
experiment was carried out using the same Data General Nova digital computer and 
driving program used for all other experiment functions. “SPOILER.FR”, the main 
program for this experiment, calls a digital I/O driver program (“DIO24”) to read 
words as fast as possible from two cards (these cards receive the digital output from 
the two Macrodyne channels, one for each LDV component). As soon as 1000 
samples (per component) have been taken by the Macrodyne processor and read by 
D1IO24, the binary output from each sample 1s converted to decimal using a decoding 
SionolumeweOr CODE LER). 

As mentioned earlier, the Macrodyne signal processor used in this experiment 
utilizes a 1 GHz counter. The number of “counts” required for the passage of 8 cycles 


is sent to the digital I/O in the form of a 10-bit mantissa and 4-bit exponent, namely 
eX 
DD De XE 


For example, if 400 “counts” were required for a particular Doppler burst, the 


output would be 


0000110010 OOI11 
Mantissa EX 


which, upon decoding, reduces to 


50 x 23 = 400 


Zi 


Since the clock counts are produced at 4 GHz (the | GHz clock is divided into 4 


phases), each clock count is 0.25 nanoseconds. Thus, 


a eX 
ig = Disp DED. 2 Pax 


is the period of the signal when 8 cycles are measured. 


Converting this to signal frequency, f., yields 


2>-EXP x 1000 


Mantissa 


Since the signal frequency, f,, 1s composed of the Doppler frequency, the shift 


frequency ‘and the mixing frequency, the Doppler frequency can be obtained from 


fp = fs 7 Tonite + fix 


Once the Doppler frequency of a particle is Known, its velocity can be determined 


from 


velocity = const X fp 


where “const” 1s determined based on the fringe spacing or direct calibration. 

With the 1000 samples from each component converted to decimal, the data can 
be reduced to average and rms velocities using the formulas found in Table [I. This 
data (both raw and reduced) is then stored for future analysis. 

While most data analysis is performed after completion of the experiment, several 
on-line analysis routines are available for use as real time performance checks. For 
example, “HISTOGRAM.FR” is a program available to first read a given number of 


samples from one of the Macrodyne channels and plot frequency versus the number of 
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samples at each frequency in histogram fashion. With this program it is easy to 
determine whether or not the data 1s conforming to the Normal Distribution (to which 
ideal data should conform). Another on-line analysis tool which ts built into the main 
program “SPOILER.FR” ts the subprogram “PLOTP.FR”. This routine gives the 
experimenter the option to have any of five velocity profiles (ic, U,)./Uo, 
VivelVoo, etc.) plotted on the screen immediately after completing an experimental 
run. This subprogram was used extensively to immediately inspect the quality of each 


run. 
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C. REYNOLDS SHEAR STRESSES AND TURBULENT STRESSES 

Table II shows the equations used to determine mean velocities, turbulent 
stresses and Reynolds shear stress throughout this experiment. The value for N was 
generally 1000 since this was the number of velocity samples taken at each 


measurement location (for each component). 
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IV. RESULTS AND DISCUSSION 


A. INTRODUCTION 

One of the purposes of this experiment was to demonstrate the benefits of the 
two-component Laser Doppler Velocimeter by emploving this twpe of svstem to 
measure the flow field of an airfoil with a deflected spoiler. In order to present the 
results of these measurements effectively, a graphics program called “PLOT-3D” was 
used to plot the data obtained. “PLOT-3D” 1s available on the VAX system at the 
Ree mes KCSearen Center, 

A data interpolation routine written by Chyang Sheng Lee of Stanford University 
was also used to interpolate the measured data into a grid svstem so that it could be 
plotted. All of the contour plots presented use the interpolation routine to interpolate 
in both the x and’ wodirections: Howewer, some of the velocity vector plots use 
interpolation in the wv direction only. This was done because, in some cases. linear 
interpolation of the data gives a false impression of what the flow field 1s actually 
doing. 


ay 
« 


In all cases, the data presented in this section is fora = 0° and6 = 30 


B. FLOW FIELD ANALYSIS 

Welocity measurements were taken at a variety of points along the upper surface 
ahead of and behind the spoiler. Figure 4.1 shows a plot of the data obtained over the 
rear portion of the airfoil. Note the low-velocity, recirculating flow region directly 
behind the spoiler. 

While the data for Figure 4.1 was interpolated in the v direction only, the data 
for Figure 4.2 was interpolated in both the x and v directions. [t can be seen in Figure 
4.2 that interpolation in the x direction vields some fairly unrealistic velocity vectors 
due to the relatively low number of measurement locations as compared to ine number 
of grid points. For example, those directly ahead of the spoiler. On the other hand, 
complete interpolation (in both the x and v directions) does serve to paint a good 
overall picture of what the flow field is doing based on a finite number of velocity 


measurements. 


y/C 


O.25 


15 


0.05 


=0,.05 


Figure 4.1 
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Flow field measured over rear portion of airfoil 
(interpolation in y direction only). 
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Figure 4.2 Flow ficld measured over rear portion of airfoil 
(interpolation in both directions). 


C. WAKE 

The wake behind the spoiler was measured with the two-component LDV system 
at numerous locations to determine the mean and rms profiles. The results from the 
wake are presented in this section. 

1. Wake Profiles 

The mean longitudinal and transverse wake profiles measured by the LDV for 
a = 0° and 6 = 30° are presented in Figures 4.3 and 4.4. The symbols are the 
measured data and the lines represent the best fit curves of the data. As expected, the 
width of the wake, the velocity defect and transverse velocity all decrease in the 
downstream direction. 

The Reynolds shear stress, longitudinal and transverse turbulent stresses 
measured by the LDV are presented in Figures 4.5 to 4.7. Obvious in all three of these 
figures is the fact that the fluctuating region appears to grow from the x/C=1.0 
(trailing edge) location to x/C= 1.25 then decrease as it approaches x/C= 1.5. The “S” 
shape profile in the Reynolds shear stress 1s typical for a two-dimensional wake 
[Ref. I:p. 29]. The longitudinal and transverse turbulent stresses show two peaks at the 
locations where the two free shear layers (corresponding to the airfoil and the spoiler) 
oscillate in the near wake. 

Figures 4.8 and 4.9 show the velocity vectors in the wake. A scaled profile of 
the airfoil is superimposed in each figure for reference. The data in Figure 4.8 was 
interpolated in the y direction only. Figure 4.9 was interpolated in both directions and 
depicts the mean wake velocity profile quite well. 

2. Contour Plots 

Figures 4.10 to 4.16 depict “contour” plots of various quantities in the airfoil 
wake. Figure 4.10 shows regions of equivalent U-magnitude (U 1s referenced to free 
stream velocity). Note that there are regions where the ratio U/Ug,) 1s greater than 
unity. This 1s due to tunnel blockage effects which cause the flow to accelerate more 
than normal around the airfoil. 

Figure 4.11 depicts the regions of equivalent V-magnitude (also referenced to 
free stream velocity). A minus sign on a contour line in Figure 4.11 implies that the V 
component of velocity is negative (downward) in that region. 

In Figure 4.12, the total velocity magnitudes in the wake are contoured. 


Again, contour lines with values greater than one arise due to blockage effects. 


34 


e3ed AdaTo 


G T=0/X 


Ge 1=0/X 





O°1=0/X 


O/A 


Figure 4.3 Longitudinal velocity profiles in the wake. 
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Figure 4.4 Transverse velocity profiles in the wake. 
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Figure 4.5 Reynolds shear stress in the wake. 
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Figure 4.6 Longitudinal turbulent stress in the wake. 
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Figure 4.7. Transverse turbulent stress in the wake. 
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Figure 4.8 Velocity vectors in the wake 
(interpolation in y dircction only). 
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Figure 4.10 Contour plot of U-magnitude in the wake. 
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Figure 4.11 Contour plot of V-magnitude in the wake. 
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Contour plot of total velocity magnitudes 
in the wake. 
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Figure 4.13 Contour plot of flow angles in the wake. 
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Figure 4.15 Contour plot of transverse turbulent stresses 
in the wake. 
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Figure 4.16 Contour plot of Reynolds shear stresses 
in the wake. 
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Figure 4.13 shows how rapidly the flow angles change in the near wake (the 
closer the contour lines, the steeper the gradient). 

Longitudinal and transverse turbulence and Reynolds shear stress are 
contoured in Figures 4.14 to 4.16. Figure 4.14 clearly shows the two regions where the 
turbulence intensity is the highest. These two regions are probably caused by 
counter-rotating vortices originating at the spoiler trailing edge (clockwise vortex) and 


the airfoil trailing edge (counterclockwise vortex). 


D. COMPARISON WITH HOT-WIRE ANEMOMETER DATA 

One purpose of this thesis was to compare the data obtained in this experiment 
with the hot-wire data obtained under similar conditions by Lee [Ref. 5]. Figures 4.17 
to 4.21 are LDV/hot-wire comparisons of longitudinal and transverse velocity, 
longitudinal and transverse turbulent stress, and Reynolds shear stress at three 
locations in the wake. 

It 1s important to mention that, in the previous experiment, no hot-wire data was 
obtained in the region of reverse flow between x/C=1.0 and x/C=1.25. Hot-wire 
anemometry 1s virtually unusable in this region due to the fact that hot-wires rectify the 
signals, which include the instantaneous forward velocity as well as the reverse velocity. 
Consequently, velocities measured by hot-wires in reverse flow regions appear higher 
than they actually are. 

Except for minor grid alignment errors in the y direction, most of the curves are 
in very good agreement. As expected however, the longitudinal turbulence stress 
curves at x/C=1.25 (Figure 4.19) and the transverse turbulence stress curves at 
x,C=1.25 (Figure 4.20) show considerable disagreement. In the longitudinal 
turbulence case, the hot-wire data shows more turbulence in the middle of the wake 
while in the transverse turbulence case, the hot-wire data shows /ess turbulence in the 
middle of the wake. This is again due to the fact that hot-wire measurements give a 
misleading velocity profile close to the trailing edge because of the reverse flow in that 
region [Ref. I:p. 70]. 

The hot-wire anemometer is quite accurate in regions where there is no reverse 
flow. However, in order to obtain reliable data on the mean velocities and turbulent 
stresses in the reverse flow regions and the near wake of this airfoil-spoiler system, 


Laser Doppler Velocimetry yields better results. 
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Figure 4.17 Comparison of LDV and hot-wire data for 
longitudinal velocity. 


50 


x/G=1 5 x/C=2.0 


x/C=1-25 





O/4 


Figure 4.18 Comparison of LIV and hot-wire data for 
transVcEse velocity. 
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Figure 4.19 Comparison of LDV and hot-wire data for 
longitudinal turbulence stress. 
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Figure 4.20 Comparison of LDV and hot-wire data for 
transverse turbulence stress. 
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Figure 4.21 Comparison of LDV and hot-wire data for 
Keynolds shear stress. 
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E. ERROR ANALYSIS 

In this experiment, 1000 velocity samples were taken by the LDV in each 
direction (LU and V) at every measurement location. The purpose in taking so many 
samples was to provide enough values in the fluctuating regions to ensure that an 
accurate statistical mean could be computed. Consequently, in the data reduction 
process all 1000 values were averaged to obtain values for U,,,, and V.,, at each 
location. The equations in Table II were then used to obtain Reynolds shear stress, 
longitudinal and transverse turbulent stress at these locations (still using all 1000 
samples). 

To determine its effect, a different data processing method was tested on the raw 
data from several runs. In this procedure, only the samples which fell within one 
standard deviation of the mean were considered (approximately 67% of the 1000 
samples) in commenting Uo; Voves Uw, % 
were employed, just as before. This new process had little effect on U 


, and V_.,.,- “gain the equauons in Table I 


ave OUt Vave 
changed considerably (from -1.3 to -11.0 at one location tested). However, the V,,.. 
values did not always decrease (one location changed from 8.7 to 16.9). Therefore, the 
changes were not due to the removal of a bias error. The longitudinal and transverse 
turbulence stress values decreased in all cases. This was expected since the biggest 
contributors to these values are the velocity terms which differ greatly from the mean. 
This new data processing procedure did not improve the data quality by 
removing any biases (based on the fact that some values of U,,, and V,,., increased 
while others decreased on using the new process). As a matter of fact, this process had 
a negative effect on the data quality in that it neglects velocity extremes which, in turn, 
lowers the apparent turbulence intensity of the flow field. For this reason, the original 


data processing procedure (consider all 1000 samples) was used in this experiment. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

A two-component LDV system was set up for the purpose of measuring mean 
and fluctuating quantities of an airfoil with a deflected spoiler. The results of this 
experiment have shown that this system provides an efficient and effective means of 
flow field measurement, including regions of reverse flow and large velocity 
fluctuations. Comparison with hot-wire data demonstrates the LDV’s superiority over 
the hot-wire in measuring in the near wake region. 

The results also clearly show two recirculating regions in the wake of the airfoil. 
One vortex originates at the trailing edge of the airfoil itself while the other originates 
at the trailing edge of the spoiler. The turbulence region is concentrated in the near 


wake and dissipates rapidly in the downstream direction. 


B. RECOMMENDATIONS 

In order to establish a more complete data base, it 1s recommended that more 
LDV measurements be obtained at other than @ = 0° and 6 = 30°. These results are 
important to the modeling of this flow field. Also, to remove the possibility of tunnel 
blockage errors, it 1s necessary to perform further studies in a large scale wind tunnel. 

When using an LDV to measure a turbulent flow field such as in this experiment, 
a large number of samples should be taken at each location. Furthermore, unless 
attempting to remove a velocity bias all of the samples taken should be considered in 


computing Uae: es etc. 
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APPENDIX 


COMPUTER PROGRAMS 


- (2) TO ADD A.M. OR P.II. 


TO THE 


SUBROUTINE STDIM (IDATE,ITIME,IVAL) 
DINENSION IDATE(S) 
DIMENSION ITINECS) 
CALL DATE (IDATE, TERROR) 
CALL TINE CITINE, PERROR) 


IF (ITIME(1).GE.13) GO TO 
LP eGimitiEtl) .&£@.12) CO 10 
CGO TO 20 

IF (ITIME(2) .CT.00) GO TO 
IF (ITIME(3).Cr.00) CO TO 
IVAL = "AM" 

CO TO 40 

IVAL = “PM" 

IF CITIME(1).LT.13) GO TO 
LEITHEC)T) = ITIMEC)) = 12 

GO TO 50 

IF (ITIMEC(1).CE.@1) GO TO 
ITIMNEC1) = ITIMEC1) + 12 

CONTINUE 

RETURN 

END 


30 
19 


30 
38 


40 


08 


217) 


"ADJUSTED" 


(1) TO CONVERT A 24 HOUR CLOCK TO A 12 HOUR CLOCK. 


TIME. 


OC 


ep 


wryly, 


3 


om 


ke) 


KORO KOKORO K KK ORK OK OK ARK ROK OK aE EK 
ORO OK OKOK OR KORO ROK ROK KOR CK KOR GK ISSR aR KK OK 


*x* 
XK 
KK 
*K 
*XK 


SPOILER.FR 
MIKE FORETMAN 


eK 
TRA 
eas 
sah 


IRS 


CK OR ORR RCE RK KES OE EK OI RIK ORR GER os 
ARK SS aS SRR OE IK CROOK RRR OK AAR OR AS ASK AS SOR SKK 


LASER DATA ACQUISITION 


PROGRAM 


RLDR SPOILER STDTM DIO24 GFLIBG 


DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DINENSTON 


IFILE(4), 
ITIHE 3), 
PRAYI C1060), 
IBINICLO6), 


INITIALIZE 


20 


DO 10 IT=1,3 
IDT(1)=90 
VARCL)=10. 
CONTINUE 


ACCEPT 
IF CINUN.LE.O) 
ACCEL 


TYPE 
CoO TV 
ACCEPT 
ACCEPT 


15 


CONL=16.635 
CON2=16.65 
VARCL) =CONL 
VAR C2) =CON2 


SUBDEVICE ADDRESSES 


ISUBL=46K 
ISUB2=4 7K 


TIME PROGIAINI 


30 


CALL STDTMCIDATE,ITIME, 


DO 30 1=1,3 
IDTC1)=IDATECT ) 
IDTC1+3)=1TInECT) 
CONTINUE 

IDT (7) =1VAL 
IDT(3)=TRON 
VARC9) = FREQ 
VARCLO) =FREQ2 


TYPE "WILL THIS BE A 
ACCEPT * 


IDT (to), 
IDAVVECS) 


"NUMBER OF SAHPLES MUST 


"MIXING FREQUENCY #1: 
"HAITNG PREQUENCY #2: 


VAR(10) 


IRAY2¢€1000) 


IBINE C16), 
V(39,20) ,X(03vu) 


-¥ C309 


"RUN NUNBER CZERO TO QUIT): 
GO TO $99 

"HUFBER OF SANPLES 
IF CNUMNM.LE.1G@O0.AND.NUMN.CE.Q) 


IVAL) 


VELI (10080) , VEL2 (1000) 


", FRUN 


PER COMPONENT: * ,NUM 
GO TO 20 
BE BETWEEN I AND 1000" 


",FREQI 
FREQ! 


(1) WAKE PROFILE" 


(2) 


BOUNDARY LAYER PROFILE 


IFCYPROF.FQ.1.OR.1TIPROF.FQ.2) GO TO 40 


TYPE "NUST ENTER IT On 


aE 
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", I PROF 


GO TO 35 
40 IFCIPROF.EQ.2) GO TO 41 
ACCEPT “X-LOCATION FOR THIS RUN CXT/C): *,VAR(3) 
GO Tv 42 
41 ACCEPT "“X-LOCATION FOR THIS RUN (X7C): *,VAR(C3) 
42 ACCEPT "*NUNBER OF Y-LOCATIONS: ",NIHAX 
IFCIPROF.EQ.2) CO TO 45 
ACCEPT"STANDARD WAKE PRUFILE Y-LOCATIONS? (1=YES,0@=NO) *°,LL 
IF(LL.EQ.1) GO TO 44 
TYPE"INPUT Y-LOCATIONS IN INCIES (LOWEST POINT FIRST):"* 
DO 43 1=1,NMAX 
WRITE (10,49) I 
ACGEPT"” ",Y(I) 
43 CONTINUE 


GO TO 30 
44 Y(1)=-1.0 
Y(2)=-0.8 
Y(3)=-0.6 
Y(4)=-0.46 
Y(3)=-0.34 
Y(6)=-@0.2 
Y¥(¢) ==0 .1 
Y(8)=0.0 
Y(9)=0.1 
Y Cro) =.2 
YCI1)=.3 
Y(12)=.4 
YC1I3)=.5 
Y(14)=.6 
YCho=.¢ 
Y¥YCre)=.6 
YOI¢)-.9 
Y(18)=1.90 
YOR) =1 ..2 
Y(20)=1.4 
GO TO 30 
45 TYPE"STANDARD BOUNDARY LAYER PROFILE Y-LOCATIONS? ° 
ACeErT" CI=YES, G-NO): * LK 


IF CLK.EQ.1)GO TO 47 
TYPL"INPUT Y-LOCATIONS IN INCHES C(CLOWEST POINT FIRST):" 
DO 46 1=1,NMAX 
WRITE(10,49) I 
AGGERT® ©. yxy(1) 

46 CONTINUE 
CO TO 390 

47 YC1)=. 
1 Gp =. 
Y(3)=. 
Y(4)=. 
Y(S)=. 
Y(6)=. 
YCGa)=. 
Y(B)=. 
Y(9)¢. 
Y(10)=1.6 
YC11)=1.1 
Yer2g)=1.2 
Y(13)=1.3 
Y(14)=1.4 
YCIS)=1.5 

49 FORMET(6X, "YC" ,12,"):") 

30 WRITECIO,129) 

125 PONE C- =. *BECIN (SCAN] .FR) »>*,2Z) 

CALL FSWAP("SCANI .SV") 


=e DONA ALON = 
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ACCEPT “DO YOU WANT RAW DATA STORED? (1=YES, @=NO) ",ITRAW 


IDT(9)=TRAW 
IDT(10)=1 PROF 


Cc 
C BEGIN TAKING DATA 


C 


N=0 


57 CALL DIOII(CISUBI , TRAY! ,1TRAY2,NUM, TER) 


C 
C DECODE 


Cc 


wo) 


we 


‘@) 


faicCan ge 
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65 


70 


735 


90 


AVE1=0. 

AVE2=0. 

DO 90 1=1,NUM 
DO 65 J=1,16 
IBINI(CJ)=1ANDCIRAY1 (1) ,2**(J-1)) 
IBIN2Z(J)=1TANDCIRAY2¢01) ,2**(J-1)) 
IFCIBINIGJ).NE.O) IBINI(J)=1 
IFCIBIN2(J).NE.0) IBIN2(J)=1 
CONTINUE 

IBN1=0 

IBN2=0 
DO 70 J=0,9 
IBN1I=IBN1+IBIN1I(J+1)*2**J 
IBNZ2=1IBN2+1BIN2(J+1)*2*%*J 
CONTINUE 

IEXP1t=0 

ITEXP2=0 
DO 75 J=10,13 
LEXP1=1EXP1+IBIN1I (J+1)*2**(J-10) 
IEXP2=1TEXP2+1BIN2(0J3+1)*2**(J-10) 
CONTINUE 


VELI (1) =CONI*(2.** (S-TEXP1 )*1000./FLOAT(CIBN1 )-40.+FREQ1 ) 


VEL2 (1 )=CON2* (2 .** (O-TEXP2)*1000./FLOAT (I BN2)-40.+FREQ2) 
AVEI=AVEI+VELI (1) 


AVE2=AVE2+VEL2 (1) 
CONTINUE 


AVE1=AVE1/NUM 
AVE2=AVE2/NUM 


CHECK VELOCITIES ARE IN LIMITS 


95 


98 


NN=N+1 
WRITE(10,92)NN,AVE2,AVE1 


FORMAT CIX, "LOCATION: “,12,4K%,"UCAVE)= “,F9.5,4X, "VCAVE)= 
IF(AVE1.GT.35.) GO TO 95 

IFCAVE2.GT.140.) GO TO 95 

GO TO 98 


WRITE FREE(10) "BAD DATA: LOCATION ",NN 


ACCEPT"RETRY OR SKIP TO NEXT LOCATION? (1=RETRY, O=SKIP): 


IF(ITRY.EQ.0)G0 TO 98 
CALL DIORM 
GO TY 37 


IFC IRAW.NE.1) GO TO 100 


STORE RAW DATA 


CALL FOPEN(S, "COUNTER" ) 
READ BINARY(S) NJ 


60 


= eA Dae aes 


", ITRY 


Gihicdsc 


108 
99 


110 


CALL FCLOSE(5) 

NJ=NJ+I 

CALL FOPEN(5, "COUNTER" ) 
WRITE BINARY(5) NJ 
CALL FCLOSE(9) 

NF ILE=NJ 

IF(N.EQ@.0) JFILE=NFILE 
IR=MODCNFILE, 10) 
IS=NIILE/100 
IL=(NFILE-IS*100)/16 
IFILEC(1I)=* RA’ 
IFILE(2)=* WD’ 

IFILE(3)= (1S+60K)*400K+(1L+60K) 
IFILE(4)= (CI R+60K)*400K 


CALL FOPEN(S, IFILE) 
WRITE BINARY(S) CIDT CI 
WRITE BINARY(S) CVARCI 
WRITE BINARY(S) CIRAY1 
WRITE BINARY(S) CTRAY2 
CALL FCLOSE(3S) 


1,10) 
1,10),NUM,YCNN) 

»1=1,NUM) 
I=1,NUM) 


ea Ae ww 
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WRITE (10,99) 

FORMAT(" ","*x*k BEGIN RMS CALCULATIONS >>",Z) 
UPNS=80. 

VRIHS=0. 

UVRUIS=0. 

DO 1!0 I=1,NUM 

URNS =URNS+ ( (VEL2 (1) )**2-(AVE2)*x2) 
VRNS=VRMS+ (CCVEL1 CI) **2-(CAVEIL )**2) ) 
UVRIiS =UVRMS+CVELI CI) *VEL2C1))-CAVEI*AVE2) 
CONTINUE 

Uhis= CURNS/FLOAT (NUM) ) 

VERMS= (VRNS/FLOATCNUMN) ) 

UVRILS= CUVRNS/FLOAT (NUM) ) 

WRITE (10,111) 

FORNHATC"+" ,"2> END RMS CALCULATIONS" ) 
N=N+1 

VC1,N)=AVE2 

V(2,N)=AVEI 

V(3,N)=URS 

V(4,N)=VRIUS 

VCS ,N)=UVERGIS 

IF(N.GE.NMAX) GO TO 120 


MOVE TRAVERSE 


CALL FSWAP(C"TRELAY.SV" ) 


ACCEPT" TRAVERSE MOVED...ANY NUMBER TO CONTINUE. 


CALL DIORI 
GO TO 37 


OBTAIN TEMNP, STAG PRES., MACII +, ETC 


120 IFCIRUN.EQ@.1609) CO TO 127? 


CALL FSWAPC"TRELAY.SV" ) 
VRITEC10,126) 


126 FORMAT C"TRAVERSE RESET TO INITIAL POSITION") 


are 


CALL FOPEN(S, "SCANI") 


READ BINARY(S) VARC4) ,VARCS) ,VARC6) ,VARCZ7) ,VAR(8) 


CALL FCLOSE(S) 


WRITE DATA INTO PERMANENT DATA FILE 
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me LAA 


999 


600 


700 


9795 


980 


CALL FOPEN(S, "PERNCOUNT" ) 

READ BINARY(S) NF 

CALL FCLOSE(CS) 

WRITE FREECIO) "LAST DATA FILE USED: ° ,MF 
WRITE FREE(1IO) “RUN NUMBER: ",TRUN 

ACCEPT "*DATA’ FILE NUMBER (NEG * TO SKIP): °,MFILE 
IF(NFILE.LT.0) GO TO 399 

CALL FOPEN(S, "PEIUICOUNT" ) 

WRITE BINARYCS) NFILE 

CALL FCLOSE(CS) 

TR=NODCHFILE, 19) 

ILety il Le7io 

IFT LEC1)=" DA’ 

IFILE(2)=° TA’ 
IFILE(3)= (CI L+601) *400K+ CIR+60K) 
IFILE(4)=6 


CALL FOPEN(S, IFITLE) 


WRITE BINARY(S) CIDT CI), 121,10) 

WRITE BINARY(S) CVARCI),1=1,10),NUM,NMAX,JFILE,NFILE 
WRITE BIRARYOS) CCVCT,J),J=2=1,NNAX),7P=1,5) 

WRITE BINARY(CS) CYCT),1T=1,NMAX) 

CALL FCLOSE(9) 

WRITE FREECIO) “ANBIENT TEMP: ",VAR(4),"° DEGF" 
WRITE FREECTIO) "ATHOSPHERTC PRESSURE: ",VARCS)," PSI" 
WRITE FREECIO) "“SPACNATION PRESSURE: " VAn(6)," Pot 
WRITE FREE C10) “FREE STREAM MACH NUMBER: ",VARC?) 

WRITE FREE (C10) "FREE STREAM VELOCITY: ",VAR(8)," M/S" 


DO 650 T=1,NrAx 

YCI)=YC1)76. 

CONTINOE 

ACCE?PT"DO YOU DESIRE A CURVE? (1=YES, O=NO): ",TCURVE 
IFCICURVE.LQ@.0) CO TO 975 

VRITE FREECLO) 

TYPE" WHAT PROFILE DO YOU DESIRE?" 


TYPE" (1) UCAVE)/UCINFINITY)" 

TYTE: (2) VOAVE)/ZUCINEINITY)" 

TYPE? (3) UCiISI7 CUCINFINITY)*«*2)" 

TYRES (4) VORMSI7CUCTNEFINE PY) #«*2) " 

ACCEPT" (9) UVGQUIS) ~<CUCINFINITY)**2) ",ISEL 
BIGK=06. 


DO 769 T=1,NNAX 

IFCISEL.LE.2) XCI)EFVCISEL, 1) /7VARCS8) 
IFCISEL.CT.2) XCT)EFVCISEL, 1)7(CVAR(CB) **2 ) 
IFCXC1) CT. BIGK) BIGN=KCI) 

CONTINUE 

CALL FOPENCS, "PROFILE" ) 

WRITE BINARY(S) NMNAX, IPROF,ISEL, BIGX 
WRITE BINARY(S) (401) ,1=21,NMAXK) 

VRITE BINARYCS) (CYCT),1=1,NMAX) 

CALL FCLOSE(S) 

CALL FSWAPC"PLOTP .SV" ) 

ACCEPT"ANY ¥* TO CONTINUE ",ANY 

TYPE "DO YOU DESIRE ANOTHER GRAPIT? " 
ACCEPT" (h=YES, O=NO):s ",IGR 
IFCIGR.EQ.1)C0 TO 600 

WRITE FREE C10) "ENSURE PRINTER IS ON" 
WRITE C12, 9SO) IDATE, ITIME, IVAL 


FORMAT CIO“, “DATES  * (2C12 , "7014, OX eee oe etn 


WRITE ERED GI) 

WRITE FREECI2)"RUN NUMBER: ", TRUN 

WRITE FREN C12) 

TF CIPROF.LQ.1) WRITE FREEC12) "WAKE PROFILE" 
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981 


IFCIPROF .£Q.2) WRITE FREEC12) "BOUNDARY LAYER PROFILE" 

WRITE PREECI2)* RXR-LOCATION: ",VAR(3) 

WRITE FREEC12) 

WRITE FREEC1I2)"NUMBER OF Y-LOCATIONS TAKEN: ",NMAX 

WRITE FREEC12) 

WRITE(C 12,981) 

FORMAT (2X, "LOCATION", 10X, "UCAVE)", 0X, "VCAVE)",10X, "UCRMS)", 


811X,"VCRHS)",11X,"UVCRMNS)") 


984 
9895 


Boo 


DO 939 I=1,NMAX 

WRITEC12,984) YC°1),V01,1),V(2,1),V03,1),V(4,1),V(5,1) 
FORHAT(C2X,F8.6,9X%,F9.4,8X%,F9.4,7K,F10.4,7X,F10.4,7X,F10.4) 
CONTINUE 

WRITE FREE(12) 

WRITE FREE(12)"NUMBER OF SAMPLES PER Y-LOCATION: ",NUM 
WRITE FREE C12) "NIX FREQ #1= ",FRFQi," MIX FREQ #*«2= ",FREQ2 
WRITE FREEC12)"VEL CONST #1= ",CONI," VEL CONST #2= ",CON2 
WRITE FREE(12) 

IFCIRAW.EQ@.@) WRITE FREE(12)"RAW DATA HAS NOT BEEN STORED" 
IFC IRAW.EQ@.1) WRITE FREE( 12) "RAW DATA IS STORED IN *RAWD’ FILES 


&JFILE," THRU ",NFILE 


WRITE FREE(12)"PERMANENT DATA FILE: DATA" ,MFILE 
WRITE FREE(12) 


WRITE FREE(12) "AMBIENT TEMP= ",VAR(4)," DEGF" 
WRITE FREEC12) "ATMOSPHERIC PRESSURE= ",VAR(S)," PSI" 
WRITE FREEC12) "STAGNATION PRESSURE= ",VAR(6)," PSI" 
WRITE FREEQ12) "FREE STREAM MACH ¥= ", VAR?) 


WRITE FREEC12) "FREE STREATT VELOCITY= ",VAR(8)," M/S" 
STOP 
END 
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*X *K 
*:K "SCANI .FR" *XK 
*KX *X 
** MIKE FOREMAN *X 
*XK *XK 


2K KK KK EK HK SK AS KK OK OK OK OEE OK OK OK OK OK KK 
2K KK KR KK KK ROK EK KK KOK RE OK OK OK KOK KK 


"SCANI.FR" TAKES TIE DATA FROM THE 6 


SCANIVALVES AND WRITES THIS DATA CIN COUNTS ) 


ON THE SCREEN AND INTO A DATA FILE. 


DIMENSION IRAY( 1200), Y(6,0:23) 


THE SCANIVALVES ARE CONNECTED TO CIIANNELS 1 THROUGH 6: 


NSTART l 
LIMIT 6 
NUN = LIMIT-NSTART+1 


ZEROIZE THE MATRIX 


DO 180 IT = 1,6 
DO 190 J = 8,23 
YCl,J> = @ 


10 CONTINUE 


SCANIVALVE RELAY INFORMATION 
IDEVICE = DEVICE NUNBER (25K FOR SCANIVALVE) 
RELAY NUMBERS: HOME = RELAY +3 
STEP = RELAY a2 
TO CONVERT RELAY NUNBER TO BINARY WORD USE: 


N=24:4(16-N) 


IDEVICE = 25K 


Nieto ;STEP = RELAY+?2 


ISTEP = 2**(16-N1) 


N2=3 ;sHOME = RELAY?3 


TWiONE = 2*x(16-N2) 
RESET ALL BITS TO ZERO 


CALL DOUTWCIDEVICE,9) 


SETUP THE A7D DEVICE WITH THE STARTING AND ENDING CHANNELS 
CALL FADSTCNSTART, LIMIT, IER) 


RESET THE SCANIVALVE TO THE HOME POSITION 
(ENERGIZE THE RESPECTIVE RELAY, THEN DEENERGIZE [T) 


CALL DOUTWCIDEVICE, IHONE) 
CALL FDELY(9) 
CALL DOUTWCIDEVICE,®) 


TAKE 200 SAMPLES AT EACII PORT 
DO 30 L=0,23 


DO 20 N=1,200 


CALL FADCTCIRAY CCN-1 )*NUM+1),NUM, TER) 
20 CONTINUE 
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ra 


Gu 


Gi wake 


cif 


C 
C 


w) 


STEP TIE SCANIVALVES TO THE NEXT PORT 


CALL NOUTW(IDEVICE,ISTEP) 
CALL FDELY(S) 
CALL DOUTWC(IDEVICE,0) 


COMPUTE THE TOTAL COUNTS FOR EACH PORT (OF EACH SCANIVALVE) 
AND STORE EACI] TOTAL IN THE ARRAY, Y. 


DO 39 N=NSTART, LIMIT 

N1=N-nNSTART+1 

DO 30 N=1,206 

YCM1,L)=YCM1 ,L)+TRAY( (N-1 ) *NUM+ (M-NSTART+#1 ) ) 
30 CONTILOE 


RETURN TIE SCANIVALVES TO THE NOME POSITION 
CALL DOUTWCIDEVICE, INONL) 


CAkL RDELYCS) 
CALL DOUTWCIDEVICE,9®) 


DISPLAY THE AVERAGE SCANIVALVE VALUES (IN COUNTS) ON THE SCREEN 
NOTE: THE AVERACE IS CONPUTED INSIDE THE WRITE STATEMENT ITSELF. 


WRITE FREEC 10)" BEGIN TEMP .SV..." 
CALE FSuart “TEMP .SV*) 


WRITE FREECIO)" EYm> TEMY.Sv >>> BEGIN’ CALPRES.SV...* 


CALL. PSMeP (*CABPRBS.SV" ) 


WRITE FREEC16)”" END CALPRES.SV >>> BEGIN NPARO.SV...” 


GALL FSWAP(*NPARO.SV") 


WRITE FROEECIO)" END NPARO.SV >>> CONTINUE SCANI.SV..." 


CALL "FOP ENS, "LEMP ” ) 
READ BINARY) TENP 
CALL FCLOSE(S) 

C&EL FOPEN(S , M@NLERES” ) 
READ BIBARY(o) CAL 
CABL FCLOSE(3) 

CAEL TOVPEN(S, “"ATMPRES®* ) 
READ BIBARY (9) ATH 
CYEL BCLOSi (5) 


ATIH=ATMN-3977768995 ;COAVERT. TO PSI 
WRITE PRLEC10) "CAL®” ,CAL 
TEHPR=TERP F499 . 6 ;CONVERT TO RANKINE 


TuTPRES= (CYC1,2)-YC1,0) )*CALZ(CYC1,1)-YC(1.0) )+ATM 
STATPRES= (Y(4 ,2) -Y(4,0) CALS CY (4. 1)-YC4,.0))+ATM 
FriACl=SQRT CABS (3.* CC CSTATPRES/TOTPRES ) ¥* (2£2.77.))71.))) 
TRY1=5.* (CC STATPRES’TOTPRES )**(-2.77.))-I1.) 
TRYZ=ABS(TRYI ) 
IFC(TRY1.EQ.TRY2) GO TO 49 
WRITE FREE(10)"STATPRES IS GREATER THAN TOTPRES !!" 
40 T1=TEHPR’ (1+ (FMACH¥*2.)75.) 
A=SQRT( 1 .4*1716.*T1)*.39048 
FVEL=A*FMACH 
GALL FOPEN(S, "SCART *) 
WRITE BIRSRY(S) TEMP, ATM, TFOTPRES, FPRMACH, FVEL 
CALL FCLOSE(S) 
REMOVE THE DEVICES 


CALL DIORM 
CALL FADRM 


SGP. 
EnD 
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"TEMP.FR" CONPUTES THE CURRENT ATMOSPHERIC 
TEMPERATURE AT THE INLET TO THE INDRAFT 
TUNNEL. THIS VALUE IS OUTPUT TO THE DISPLAY 
SCREEN AND ALSO WRITTEN INTO A FILE CALLED 
was Wa Dw | eer 
DIMENSION IRAYC200) 
THE TENPERATURE PROBE INTERFACES WiTH CHANNEL 8 
NSTART = 8 
LINIT = 8 
NOM = LIMNIT-NSTART+ I 
SET UP THE DEVICE 
CALL FADSTCNSTART, LIMIT, 1 ER) 
TAKE 200 SAMPLES 
DO 10 I[=1,200 


CALL FADCTCIRAYCIT) ,NUN, IER) 
10 CONTINUE 


COMPUTE THE TOTAL NUMBER OF COUNTS FOR THE 200 SAMPLES 
TEHP=0. 
DO 20 1=1,200 
JEnP=TENP+ITRAY CT) 

20 CONTINUE 

CALCULATE TITE AVERAGE OF THE 200 SAMPLES 
TENMP=TENP/200. 


THE TEMPERATURE SIGNAL ENTERS THE AMPLIFIER IN VOLTS. 

THE ANPLIFIER GIVES THE SIGNAL A 6ODB GAIN (FACTOR 

OF 1000). THEREFORE, THE SIGNAL GOES TO TIE A’“D CONVERTER WITH DIMENSIONS 
OF MNILLIVOLTS. AFTER CONVERSION, Tike FACTOR 1638.3 COUNTS/MILLI VOLT 

IS NECESSARY TO CONVERT THE COMPUTER UNITS (COUNTS) BACK TO MILLIVOLTS: 


TEMNP=TEHP/1638.3 ;CONVERT COUNTS BACK TO MILLIVOLTS 
THE "TEMP" HAS NOW BEEN CONVERTED TO MILLIVOLTS 
HOWEVER, THE TEMPERATURE EQUATION REQUIRES THAT 
THE "“TENP" BE IN MICROVOLTS: 


TEMP=TEMP*1000. 
BECKMAN COEFFICIENTS 


CO=3 .1922787E+01 
C1l=3.9497620E-02 
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C2=-2.4786655E-07 
CALCULATE TEMPERATURE 
TEMP =CO+C1*TEMP+C2*TEMP*«2 
WRITE TEMP TO SCREEN AND FILE "TEMP" 
WRITE FREE (10) 
WRITE FREE (10) “TEMPERATURE (DEGF)= “,TEMP 
CALL FOPEN (5,"TEMP") 
WRITE BINARY (5) TEMP 
CALL FCLOSE (5) 
DISCONNECT THE DEVICE 
CALL FADRM 


STOP 
END 
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"CALPRES.FR" THIS PROGRAM COMPUTES THE 
CALIBRATED PRESSURE TO THE SCANIVALVES. 


DIMENSION ITRAY(200) 


THE CALIBRATION PRESSURE IS CONNECTED TO CHANNEL 7. 


NSTART = ? 
LINIT = 7 
NOM = LIMIT-NSTART+1 


SET UP THE A/D DEVICE: 


CALL FADSTCNSTART, LIMIT, LER) 


TAKE 200 SAMPLES AT CITANNEL 7 


DO 19 T=1,200 


CALL FADCTCIRAYCT) ,NUM, TER) 


10 CONTINUE 


CONPUTE THE AVERAGE OF TITE 200 SAMPLES 


Y=0. 

BO 20 1=1,200 

Y=Y+FLOAT CTRAY CT) ) 
20 CONTINUE 

Y=Y/269. 


CALIBRATION PRESSUINE IN INCIES OF HG = VOLTS*10 


CAL= (CY71638.3):«10. 
CALIBRATION PRESSURE IN PSI 


CAL=CAL* (3377 .76895.) 
WRITE FREE(10) CAL 


WRITE THE CALIBRATION PRESSURE TO THE SCREEN AND TO 
A TEMPORARY FILE CALLED "CALPRES". 


WRITE FREE(10) 
WRITE €10,29) CAL 


29 FORMAT ("CALIBRATION PRESSURE: 
CALL FOPEN (5, “CALPRES" ) 


WRITE BINARY (5) CAL 
CALL FCLOSE (3) 


REMOVE THE DEVICE 
CALL FADRM 


STOP 
END 
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INCHES OF HG * (337776895) 


rol) 


C "PLOTP.FR" IS CALLED BY "SPOILER.FR* TO PLOT SELECTED VELOCITY 
C PROFILES. 


C 
DIMENSION X(25), Y(25) 
DIMENSION IST1(°3),1ST2(°8),1ST3(8) , IST4(8) 
DIMENSION ISTS(8),1STE6(9), 1IST7(21),1ST8(31) 


CALL CIRS1Z(4) 

DATA IST1789 ,47 ,677 

DATA IST2/783,40,65,66,69,41 ,47,85/ 

DATA IST3/86,40,65,86,69,41,47,85/ 

DATA IST4/683,40,82,77,83 ,41,47 ,857 

DATA IST9/786,49 ,82,77 ,83,41,47,857 

DATA IST6/895 ,86,40,82,77 ,83,41 ,47,85/ 

DATA IST?/8&7 ,63,75, 69,32 ,86,69,76,79,67,73,84,89,32, 
*689,82,79,70,73,76,697 

DATA IST8/’66,79,85,78,68,653,82,89 ,32,76,65,89,69,82,32, 
*80,69,76,79,6¢0,73,84,89 ,32,80,62,79,70,73,76,69/ 


CALL FOPEN(CS, "PROFILE" ) 

READ BINARY(S) NNAX, IPROF, ISEL, BIGX 
READ BINARY(S) (CXC1),1=1,NMAX) 

READ BINARY(OS) CYCT),1=1, NMA) 
G3LUFCcLrose cs) 

IF C(BIGX.LT..3) BIGX=.5 

WRITE FREECIO)YC1),YCNMAX) , BIGX 
ACCEPT"ANY NUMBER TO CONTINUE: ",IRAND 


CALL INITT( 960) 

GALL BINTTT 

CALL CrntslZ(3) 

CALL TERII(C3 , 1023) 

CALL SLITAX(1G60,900) 

CALL SLITIYC130,690) 

CALL AFRITCS) 

CALL YFITUICS) 

CALL XTICS(10) 

GALE YIICS C190) 

CALL DLITIX(6.6, BIGX) 

CALL DLINHYCYC1) , YCNMAX) ) 

CELL PIS CNA) 

CALL CHECACX,Y) 

CALE DSPLAY (X,Y) 

CALL CHUSTIZ(3) 

CALL TOVALS (10,525) 

CALL VLABEL (3,IST1) 

CAI.L MOVABS (300,725) 
IF(IPROF.E@.1) CALL HLABEL(21,1ST?7) 
IFC IPROF.EFQ@.2) CALL HLABEL(3I1,IST8) 
CALL MOVABS (428,75) 

IFCISEL.EQ@.1) CALL HLABEL(8, IST2) 
IFCISEL.EQ.2) CALL HLABEL(8,IST3) 
IFCISEL.EQ.3) CALL HLABEL(8, IST4) 
IFCISEL.EQ.4) CALL NLABEL(8S,ISTS) 
IFC ISEL.EQ.3) CALL NLABEL(9,1IST6) 
CALL FINITIT(0,750) 

CALL CHRSIZ(4) 

STOP 

END 
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IISTOGRAM.FR RECORDS DATA FROM THE LDV 
COUNTER AND PLOTS A HISTOGRAM OF THIS DATA. 


REQUIRED SUBROUTINES: DIO22 eFLIBG 


DINCNSION IRAY( 1000), IBINCI6) , FREQ(1000) 
DINENSION XC101),YC1OL),VELC1000) 


TYPE “INPUT SUBDEVICE CODE --- TWO DIGITS ONLY: “ 
TYPE. CHANNEL 1 = 46" 
PYRE e: CHANNEL 2 = 47" 


READ(1I1,9) ISUB 
S FORILVY (O02) 
ACCEPT "VELOCITY CONSTANT? ",CON 
ACCEVT "MIXING FREQUENCY? ", FREQM 
10 ACCEPT “NUIBER OF SAHPLES (0 TO STOP) ?" ,NUM 
Ik (NUM.CT.1600.OR.NUI.LT.G) CO TO 10 
IF (NUH.EQ.9) GO TO 999 
CALL DIOTICISUB, IRAY, NUT, PER) 


DO 100 J=1,NUN 
DD 50 IT=1,16 
IBDINCI)=LANDCIRAY( J) ,2%*CI-1)) 
00 IFCIBINGI).NE.O) IBINCI)D=1 
IBN=9 
DD 60 1=0,9 
IBN=IBNtIBINCI +1 #2] 
60 CONTINUE 


IEXP=90 
DO 70 1=10,13 
PEXP=TEXP+IBINC] +1) *2:&* (1-10) 
79 CONTINUE 


FREQ(J)=32000. “FLOAT CI BN:2:% T EXP ) 
VEL(CJ) =CON&* CFREOQ (J )-40.+FREQM) 
100 CONTINUE 


AVER=0. 
AVLE=0. 
DO 110 IT=1,NUM 
AVER=AVER+FREQCT ) 
AVE=AVE+VELCI ) 
110 CONTINUE 
AVER=AVER/NUM 
AVE=AVE/NUEI 
WRITE FREECIO) "AVERAGE DOPPLER FREQUENCY= ",AVER, "MHZ" 
WRITE FREECIO) "AVERAGE VELOCITY= ",AVE," M78" 
ACCEPT "DO YOU WISI A HISTOCHRAN? CYES=1, NO=0)",AA 
IFCAA.EQ.0) CO TO 999 
MMN=TFEXCAVE) 
JJ=NN-3S 
IFCJIJ.LT.0) JJ=@ 
KK=JJ+190 
LL=(JJ*10)-1 
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150 


999 


Z=0. 

DO 200 171,101 

K=1+LL 

X(1)=(FLOAT(K))716. 

Y(I)=0. 
DO 150 J=1,NUM 
TI=IFIX(10.*VEL(J)) 
IFCII.EQ.K) YCI)=FYCI) +1 
CONTINUE 

WRITE FREE (10) KCI),YCI) 

IW=IW+IFIXCYCT)) 

IFCYCI).GT.Z) Z=YC1) 

CONTINUE 

Z=FLOAT (CLF IXCZ/100. )*100)+100) 

IF(Z.GPF.1069.) Z=1000. 

AB=FLOAT (JJ) 

AC=FLOAT (KK) 

WRITE FREE (10) AB,AC,Z 


CALL FOPEN (5, "HISTO" ) 

WRITE BINARY (5) AB,AC,Z 

WRITE BINARY (5) CXC€1), 1=1,101) 
WRITE BINARY (5) CYCI), 1=1,101) 
CALL FCLOSE (5) 

CALL FSWAP C°HISTOGRAPII.SV" ) 
WRITE FREECIO) 

WRITE FREE(C10) 

WRITE FREECIO) 

WRITE FREECIO) 


WRITE FREE (10) "“ TOTAL SAMPLES PLOTTED: 
CALL DIORM 

STOP 

END 
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C "“HISTOGRAPH.FR" USED BY HISTOCRAM.FR TO PLOT 
C A HISTOGRAM OF LASER DATA. 
C 

DIMENSION X(101),YC1O01) 

DIMENSION ISTIC?@), IST2(°21), IST3C1S) 


O 


DATA IST1783,65,77,80,76,69,837 

DATA 1S8ST2772,73,83,84,79.71 ,82,65,77, 
*32,79,70,32,76,68,86,32,68,65,84,65/7 

DATA IST3776,82,69,81 ,85,69,78,67,89,32,40,77 ,72,90,417 


'S) 


CALL FOPEN(S, "HISTO" ) 
READ BINARY (5) AB,AC, 
READ BINARY (35) (XH), 
READ BINARY (5) CYCH), 
CALL FCLOSE (5) 


,101) 


Zz 
I=1 
T=1,101) 


@) 


CALL INITT( 960) 

CALL BINITF 

CALL TERN(3, 1023) 
CALL VBARST(C15.6,3) 
CALL SLINN(100,900) 
CALL SLINY(190,6590) 
CALL XFRIH(2) 

CALL YEFRIHCS) 

CALL XTICS(10) 

CALL YFICS(10) 

CALL AMPCS (10) 

CALL YHres <2) 

CALL DLIMNA(CAB,AC) 
CALL DLIMY(0.0,Z) 
CALL NPTS CIOL) 

CALL CHECK(X,Y) 

CALL DSPLAY(X,Y) 
CALL NOVABS (10,525) 
CALL VLABEL(7,IST1) 
CALL MOVABS (350,725) 
CALL HILABEL(21,1ST2) 
CALL MOVABS (420,75) 
CALL HWLABELCIS, ISTS3) 
CALL NUVABS (900,150) 
CALL DPRWABS (900, 650) 
CALL CHURSIZ(4) 

CALL FINIFFCO,750) 
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